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Abstract 

In a diquark-diquark-antiquark picture of pentaquarks, we use two interpolating currents to 
(^ ■ calculate the mass of the recently measured r, state in the framework of QCD sum rules. We 

o. 

CN . show that, even though yielding similar values for m=-- (and close to the experimental value), 

-)— > ' 

^ \ these currents differ from each other in what concerns the strength of the pole, convergence of the 

. -^ ■ OPE and sensitivity to the continuum threshold parameter. 
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the central question to be 



I. INTRODUCTION 

After the recent discovery of pentaquark states 

addressed now concerns the structure of these new baryons. They could be: a) uncorrelated 

quarks inside a bag p]; b) a K — N molecule bound by a van-der Waals force |l|; c) a 

n 

^^K — N" bound state in which uud and us are not separately in color singlet states 8]; d) 
a diquark-triquark (ud) — (uds) bound state [9| and e) a diquark-diquark-antiquark (DDA) 
state. This last one was advanced by Jaffe and Wilczek (JW) 'l(|] and is quite appealing 
because it can explain two unusual features of pentaquarks: their small mass and decay 
width. Instantons generate strong attractive quark-quark interactions with the formation of 
low mass diquarks, which, in turn, lead to relatively low mass pentaquarks. This was verified 
in 11 1 , [12] and |l3|. Moreover, in [l^, it was shown that the DDA configuration may lead to 
strongly suppressed transition amplitudes (to meson-baryon states) for a reasonable choice 
of its spatial wave function, namely, two separated extended diquark balls overlapping only 
partially and the antiquark at the center of the system. 
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Pentaquark configurations can be implemented in QCD sum rules (QCDSR) 
and in lattice QCD |l7|, |l8| by a proper choice of the interpolating current. A current 
for configuration c) has been proposed by Zhu J8[ and for configuration e) three different 
currents were proposed in |ll|, [l^ and 13| |. 

Given the impulse of this field and the prospects of new measurements, we may expect 
that the "pentaquark wave" will last still for a long time. A great effort will be devoted 
to understand the structure of these objects. For the QCDSR community (and for lattice 
QCD studies as well) this means that more attention will have to be given to the properties 



of the interpo 



atin g cu rrents. At this point, we can already compare the three calculations 



p..esented . U, Q a„d Q f„. tKe .as. of the 0^ A.thou,. they present difeent 
implementations of the DDA scheme, they produce very similar results for the pentaquark 
mass. This indicates that a more careful analysis has to be performed. 

The purpose of the present work is twofold. We will use QCDSR to calculate the mass 
of the recently measured H and we will also take the opportunity to perform a careful 



comparison of the results obtained with two different currents, which we wil 
Current I was introduced by us in [ll| and current II is the one proposed in 



call I and II. 
. As it will be 



seen, both currents give approximately the same mass for the S but have different sensi- 



tivity to the continuum threshold parameter, different convergence of the operator product 
expansion (OPE) and also different strength of the pole, with respect to the continuum. 
This comparative study, which will later be extended to the currents suggested in [ISf and 



also in 19[ , will give us a better understanding of these currents and may lead to the choice 



of the best current. 



II. CURRENTS AND CORRELATION FUNCTIONS 



In the QCDSR approach 
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16|, the short range perturbative QCD is extended by 
an OPE expansion of the correlators, which results in a series in powers of the squared 
momentum with Wilson coefficients. The convergence at low momentum is improved by 
using a Borel transform. The expansion involves universal quark and gluon condensates. The 
quark-based calculation of a given correlator is equated to the same correlator, calculated 
using hadronic degrees of freedom via a dispersion relation, providing sum rules from which 
a hadronic quantity can be estimated. The QCDSR calculation of hadronic masses centers 
around the two-point correlation function given by 

n(g) = i fd^xe"^''{0\Tr]{x)ri{0)\0) , (1) 

where ri{x) is an interpolating field (a current) with the quantum numbers of the hadron we 
want to study. 

A. Current I 

Following the diquark-diquark-antiquark scheme, we can write two independent interpo- 
lating fields with the quantum numbers of S : 

Vi{^) = -j=eabc{d^{x)C'j5Sb{x))[d^{x)C-f5Se{x) + sf (x)C754(x)]CmJ(x) (2) 

V2{x) = -^eaU€{x)Csk{x))[d^{x)Cse{x) + s^ {x)Cd,{x)]Cu^ (x) (3) 

where a, b, c and e are color index and C = —C^ is the charge conjugation operator. 

As in the nucleon case, where one also has two independent currents with the nucleon 
quantum numbers |20,|2l|, the most general current for S is a linear combination of the 
currents given above: 

ri{x) = [tr]i{x) +r]2ix)], (4) 



with t being an arbitrary parameter. In the case of the nucleon, the interpolating field with 
t = — 1 is known as lofi:e's current |20|. With this choice for t, this current maximizes the 
overlap with the nucleon as compared with the excited states, and minimizes the contribution 
of higher dimension condensates. In the present case it is not clear a priori which is the best 
choice for t. This will be studied in what follows. 

B. Current II 



We also employ the following interpolating field operator for the pentaquark 
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r^(x) = e''''e'''^e'f'^{s];ix)Cdt{x)}{s^ix)C^Mx)}Cu^gix) (5) 

It is easy to confirm that this operator produces a baryon with J = 1/2, / = 3/2 and 
strangeness —2. The parts, S^{x) = e"'"^s^(a;)C75(if,(a;) and P'^{x) = e"-'^'^s^{x)Cdb{x), give 
the scalar S (0"*") and the pseudoscalar P (0^) sd diquarks, respectively. They both belong 
to the anti-triplet representation of the color SU(3). The scalar diquark corresponds to the 
the I = 1/2 sd quark with zero angular momentum. It is known that a gluon exchange force 
as well as the instanton mediated force commonly used in the quark model spectroscopy 
give significant attraction between the quarks in this channel. 

C. Mass sum rules for the S 

Inserting Eq. (^ into the integrand of Eq. (^ we obtain 

{0\Tr]{x)rj{0)\0) = t'Uu{x) + t {Uu{x) + U2i{x)) + Haslx) (6) 

Calling Fi = 75 and r2 = 1 we get 

U,,{x) = {0\Tr],{x)r]^{0)\0) = eabcea'b'c'CS,J{-x)C {- Tr [riS',,,{x)r,CSl,{x)C] 

X Tr [r,s:A^)r,csJA^)C]+ Tr [r,sUx)r,csJA^)cr,s:Ax)r,csJA^)C] 

+ Tr [T,StAx)T,CSjAx)CT,SUx)T,CSl,{x)C] - Tr [T,SUx)T,CSUx)C] 

X Tr [r,s:Ax)rjCSjAx)C] - Tr [r,Sab'{x)r,cs:J,{x)cr,s,e'{x)r,cs'J,{x)C] 
- Tr [r,sUx)r,csJAx)CT,s^eA^)r,cs:A^)C] + Tr [r,sUx)r,cs:Ax)c] 

X Tr [T,S,Ax)T,CS:,Ax)C] + Tr [T,SUx)T,CSJAx)C] Tr [T,S:Ax)T,CSJAx)C]} ,(7) 



where Sab{x) and S^fj{x) are the hght and strange quark propagators respectively. 

In order to evaluate the correlation function 11 (g) at the quark level, we first need to 
determine the quark propagator in the presence of quark and gluon condensates. Keep- 
ing track of the terms linear in the quark mass and taking into account quark and gluon 
condensates, we get J22] 






where we have used the factorization approximation for the multi-quark condensates, and 



we have used the fixed-point gauge 



22|. 



Inserting (jHJ into (|7j) we obtain a set of diagrams which contribute to the OPE side of the 
correlation function. In the case of the current I, Fig. 1 shows the classes of diagrams which 
we are considering. In each diagram, except the perturbative one, all possible permutations 
of the lines are implicit. 

Lorentz covariance, parity and time reversal imply that the two-point correlation function 
in Eq. (^ has the form 

n(g)=n,(g2)+n,(g2)^. (9) 

A sum rule for each scalar invariant function Hi and Ilg, can be obtained. As in ref. J8|, 
in this work we shall use the chirality even structure ng(g^) to obtain the final results but 
we shall also discuss the chirality odd structure Hi, showing its limitations for our present 
purposes. 

The phenomenological side is described, as usual, as a sum of pole and continuum, the 
latter being approximated by the OPE spectral density. In order to suppress the condensates 
of higher dimension and at the same time reduce the influence of higher resonances we 
perform a standard Borel transform [l5|: 

Ii{M')^ lim l(Q2)n+iL M"n(Q2) (10) 

n,Q2^oo n! y ttQ / 

(Q^ = — g^) with the squared Borel mass scale M^ = Q^ /n kept fixed in the limit. 

For current II we repeat the steps mentioned above, substituting (jSJ into ([T)), making use 
of the expansion (jH}, picking up the terms multiplying the structure ^ and finally performing 



a Borel transform pU|) . In this case the diagrams that have to be considered are shown in 
Fig. 2. 

After Borel transforming each side of UqlQ"^) and transferring the continuum contribution 
to the OPE side we obtain the following sum rule at order mg: 



mi. fSn 

,2 - ' 



Aie"^=/ e-^p1{s)ds. (11) 



/o 

where i{= I, II) refers to the current employed and the spectral densities, up to order 6 are 
given by: 

s ex s 

■5^ / - 2 - 2\ ■S^ 

+ Co o n — 7 < ss > + < 00 > + C4 . „ , < ss >< qq > 
— C5— — -— - — rm„ < sOoCr.Gs > 



^'^4 ,,„,„i^ fi ^s < qgsCr-Gq > 61n(- 



'4!3!2i%6-- ^..^-". ' V"""A^^^^ 2;' 



(12) 



with ci = 5t2 + 2t + 5, C2 = (1 - tf, C3 = (t + 1)^ C4 = t^ - 1, C5 = t'^ + 22t + 1, Aqcd = HO 
MeV and 

s s s a 

P''^'^ = 5!5!2i077r8 + 5!3!2^^'"^ < '' ^ +5!3!2i07r6 < f ^' > 

+ . . „ ^ m. < sOoCr.Gs > . (13) 

4!3!2%6 '^ y^ \ i 

To extract the H mass, 771=, we take the derivative of Eq. ()14p with respect to M~^ and 
divide it by Eq. (fT^ . Repeating the same steps leading to (fT^ . (fT^ and (fT^ for the chirality 
odd structure ni(g^) we arrive at 

A|m=e~"M^= / e'M^Pj^(s)(is. (14) 

JO 



where 



and 



P^i (^) = -^1 51412107,6 < ^5 > + ^1 4,3,2127^6 < ^^-^-^9 > (15) 

P'^' ^'^ = ~5!4!2^ ^^^^ + 4!3!2%^ < ^"^^^'^^ ^ ^^^^ 

6 



In the complete theory, the mass extracted from the sum rule should be independent of 
the Borel mass M^. However, in a truncated treatment there will always be some dependence 
left. Therefore, one has to work in a region where the approximations made are supposedly 
acceptable and where the result depends only moderately on the Borel variables. 

III. EVALUATION OF THE SUM RULES AND RESULTS 

Before proceeding with the numerical analysis, a remark is in order. A comparison 
between results obtained with different currents is more meaningful when they describe 
the same physical state, i.e., those with the same quantum numbers. Concerning spin, all 
currents considered in our work have the same spin (= 1/2). Concerning the parity, the 
situation is more complicated. In QCD sum rules, when we construct the current, it has a 
definite parity. Current I has parity P = — 1 and current II has parity P = +1. However, 
currents can couple to physical states of different parities. As well discussed in J2^ , in order 
to know the parity of the state in QCDSR, we have to analyze the chiral-odd sum rule. If 
the r.h.s of this sum rule (containing the spectral density coming from QCD) is positive, 
then the parity of the corresponding physical state is the same as the parity of the current. 
If it is negative the parity is the opposite of the parity of the current. Performing this 
analysis we might determine in both cases the parity of the state. However it turns out 
that for both currents, in the chiral-odd sum rule the OPE does not have good convergence, 
i.e., terms containing higher order operators are not suppressed with respect to the lowest 
order ones. This sum rule is thus ill defined and nothing can be said about the parity of the 
state. The comparative study of the currents is still valid because we can compare other 
properties of these currents. A possible outcome of this study might be that one current has 
defects, which are so severe that we are forced to abandon it. If this turns out to be case, 
the determination of the parity of the associated states becomes irrelevant. 

Having presented the main formulas in the last section, the next step will be to introduce 
numerical values for the masses and condensates, choose reasonable values for the free (or 
partially constrained) parameters, which are the continuum threshold, the Borel mass at 
which the mass sum rule is evaluated and, in the case of current I, the value of t. As a 
result we obtain values for mH. In doing these calculations we must remember that it is not 
enough to obtain a pentaquark mass consistent with the experimental number. There is a list 



of requirements that must be fulfilled: 

i) the physical observables, such as masses and coupling constants, must be approximately 

independent of the Borel mass (this is the so called Borel stability). 

ii) the right hand side (RHS) of the sum rule (J14j) must be positive, since the left hand side 

(LHS) is manifestly positive. 

iii) the operator product expansion (OPE) must be convergent, i.e., the terms appearing in 

(fT^ and (fT^ must decrease with the increasing order of the operator. 

iv) the pole contribution must be dominant, i.e., the integral in (jT^ must be at least 50 % 

of the integral over the complete domain of invariant masses (sq -^ cxo). 

v) the threshold parameter Sq must be compatible with the energy corresponding to the first 

excitation of an usual baryon. 

vi) the current-state overlap Ah must be sizeable. The larger is Ah th better is the current. 

As it will be seen the above conditions are not easy to be simultanously satisfied and may 
be useful in discriminating between different currents. 

In the numerical analysis of the sum rules, the values used for the condensates are: 
{qq) = -(0.23)3 GeV^ {ss) = 0.8{qq), < sg^a-Gs >= ml{ss) with ml = 0.8 GeV^ and 
{qIG'^) = 0.5 GeV"^. We define the continuum threshold as: 

So = (1.86 + Af GeV^ (17) 

A. current I 

We evaluate our sum rules in the range 1.0 < M^ < 4.0 GeV^. The results are shown in 
Table I. The first column in Table I gives the values of t considered here. We notice that 
some values were excluded, as for example t = 0, because they would lead to a violation of 
condition ii) above. The simbol S stands for a current composed by scalar diquarks only, 
i.e., r]i. The motivation for studying this current is to verify if it gives a smaller mass for the 
pentaquark than those obtained with other currents. According to the instanton description 
of diquark dynamics, this should be the case. In the second and third columns of Table I 
we list the values adopted for the strange quark mass and continuum threshold respectively. 
The fourth column shows rris and the Borel mass squared at which the particle mass was 
extracted. The fifth and sixth columns show the ratios of terms containing dimension 4 
(quark and gluon condensates) and 6 (mixed condensates) operators and the perturbative 
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term respectively. Ideally these numbers should be smaller than one, the second being 
smaller than the first. Finally, the last column shows the strength of the pole. In practice 
this quantity is obtained dividing the integral in the right hand side of (J14j) by the same 
integral with the upper limit Sq = oo. This percentage must be as large as possible, but 
usually we accept values around 50 %. 

As it can be seen in Table I, it is hard to satisfy conditions i) - vi) simultaneously. In 
particular, when we have a very good OPE convergence, the strength of the pole is very weak 
and vice versa. We have to look for a compromise and we believe that the choices indicated 
in bold face are the best. For same choices for rus and A we observe that the masses found 
with scalar diquark currents, S, are only slightly smaller than the others. This means either 
that instanton dynamics was not captured by our choice of currents and diagrams or that 
the interaction between the pseudoscalar diquarks (included in the mixed currents) is more 
attractive than expected. The connection between our approach and instantons deserves 
further investigation. 

In order to further illustrate these results, we consider the parameters of the first line in 
bold face, rris = 0.10 GeV and A = 0.44 GeV, and contruct. Figures 3, 4 and 5, showing the 
Borel mass dependence of m=, of the OPE terms (in absolute value) and of the percentage 
of the pole contribution respectively. For these choices the value of the current-state overlap 
is: 

A? ~ 5.4 X 10-^ GeV^3 (18) 

In view of the variety of results presented in Table I and in the figures, it would be some- 
what artificial to invent a procedure to determine the average value of the S mass and its 
corresponding error. If we simply take the average over all values presented we arrive at: 

ms = 1.87 ±0.22 GeV (19) 

Restricting ourselves to the parameter combinations which satisfy the requirements i) - vi) 
and taking the average we obtain: 

m^ = 1.85 ±0.05 GeV (20) 

These numbers are very close to the experimental value. 



B. current II 

Using the same numerical inputs quoted in the last subsection we evaluate now the sum 
rules obtained with current II. The results are shown in Table II, which has the same format 
as Table I except for the fact that, now, there is no t parameter. The same comments made 
for the results of Table I apply here. Choosing the parameters of the first line in bold 
face, rris = 0.10 GeV and A = 0.24 GeV, we present in Figures 6, 7 and 8 the Borel mass 
dependence of m-E, of the OPE terms and of the percentage of the pole respectively. As it 
can be seen, with current II we tend to overestimate m^, unless very low threshold parmeters 
or quark masses are used. Besides, the pole contribution is always smaller than 40 %. On 
the other hand, the Borel stability seen in Fig. 6 is remarkable. This suggests that we could 
choose a lower value for the Borel mass, thereby increasing the pole contribution without 
significantly changing rris- In order to extract from Table II some value for the pentaquark 
mass, we repeat the procedure applied in the last subsection. The simple average gives: 

ms = 2.00 ±0.16 GeV (21) 

and the average over the best parameter choices leads to: 

ms = 1.88 ±0.04 GeV (22) 

These numbers are somewhat high but still compatible with data. Finally, for these param- 
eters the current-state overlap is: 

A? ^ 1.3 X 10"^ GeVi3 (23) 

Comparing (fTHj) and ()23|) we observe that the coupling of current I to the the S state is 
four times larger than the coupling of current II to this state. This speaks in favor of current 
I. 

IV. SUMMARY AND CONCLUSIONS 

In conclusion, we have presented a QCD sum rule study of the H pentaquark mass 
using a diquark-diquark-antiquark scheme. We have employed two possible currents (I and 
II) which implement this picture. Several parameter choices were considered for the currents. 

10 



We could conclude that for both currents it was possible to find a set of parameters and a 
Borel window so that conditions i) - vi) were fulfilled and a reasonable value for the mass 
was found. The ratio pole/total, presented in the tables is still small and probably is the 
most important drawback of the studied currents. 

At the present stage we may say that the pentaquark S can be reasonably well described 
by QCD sum rules and the value of its mass is given by ()20|) or by (J22I)- These numbers can be 
considered the same within the uncertainties. Concerning the currents, there are differences 
in the Borel stability, in the sensitivity to the strange quark mass, in the continuum threshold 
and in the their couplings to the S state. 
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t 


m, {GeV) 


A (GeV) 


m^{Ml) (GeV) 


Dim 4/Pert 


Dim 6/Pert 


Pole/Total 


-1 


0.15 


0.44 


1.95 (3.5) 


2.42 


3.67 


2.1% 


-1 


0.15 


0.44 


1.66 (1.0) 


0.36 


'~7 '~7 

i . i 


80% 


1 


0.15 


0.44 


2.09 (3.5) 


0.96 


0.53 


0.28 % 


1 


0.15 


0.44 


1.82 (1.0) 


1.44 


1.11 


40% 


s 


0.15 


0.44 


1.94 (3.5) 


2.42 


11.6 


4.8% 


s 


0.15 


0.44 


1.65 (1.0) 


3.62 


23.9 


87% 


-1 


0.12 


0.44 


1.95 (3.5) 


0.24 


3.8 


2.1% 


-1 


0.12 


0.44 


1.66 (1.0) 


0.36 


8.0 


80% 


1 


0.12 


0.44 


2.09 (3.5) 


0.76 


0.42 


0.32 % 


1 


0.12 


0.44 


1.86 (1.0) 


1.15 


0.89 


41% 


s 


0.12 


0.44 


1.94 (3.5) 


1.95 


11.6 


4.8% 


s 


0.12 


0.44 


1.65 (1.0) 


2.93 


24.1 


87% 


-1 


0.10 


0.44 


1.95 (3.5) 


0.24 


3.89 


2.1 % 


-1 


0.10 


0.44 


1.66 (1.0) 


0.36 


8.2 


80% 


1 


0.10 


0.44 


2.09 (3.5) 


0.64 


0.35 


0.34 % 


1 


0.10 


0.44 


1.88 (1.0) 


0.96 


0.74 


41 % 


s 


0.10 


0.44 


1.93 (3.5) 


1.64 


11.7 


4.8% 


s 


0.10 


0.44 


1.64 (1.0) 


2.47 


24.3 


87% 


1 


0.09 


0.44 


2.09 (3.5) 


0.57 


0.31 


0.36 % 


1 


0.09 


0.44 


1.90 (1.0) 


0.86 


0.67 


42 % 


1 


0.09 


0.34 


2.00 (3.5) 


0.68 


0.40 


0.23 % 


1 


0.09 


0.34 


1.81 (1.0) 


0.97 


0.78 


34% 



Table I: Results for the H obtained with current I. 



13 



rris (GeV) 


A (GeV) 


m^{Ml) (GeV) 


Dim 4/Pert 


Dim 6 /Pert 


Pole/Total (%) 


0.15 


0.24 


2.03 (3.5) 


0.02 


0.44 


0.11 % 


0.15 


0.44 


2.16 (3.5) 


0.02 


0.26 


0.35 % 


0.12 


0.14 


1.88 (3.5) 


0.35 


0.46 


0.10 % 


0.12 


0.14 


1.93 (1.0)2 


0.46 


0.76 


16% 


0.10 


0.24 


1.93 (3.5) 


0.47 


0.29 


0.22 % 


0.10 


0.24 


1.87 (1.0)2 


0.64 


0.52 


28.7 % 


0.09 


0.24 


1.84 (1.0)2 


0.76 


0.46 


31.1 % 


0.09 


0.34 


1.91 (1.0)2 


0.67 


0.39 


39.0 % 



Table II: Results for the 



obtained with current II. 
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FIG. 1: Classes of diagrams which contribute to the OPE of current I up to dimension 6: a) 
perturbative; b)-d) quark condensate; e) gluon condensate; f)-h) 4-quark condensate; i)-o) mixed 
condensate. 
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FIG. 2: Diagrams which contribute to the OPE of current II: a) perturbative; b) quark condensate; 
c) gluon condensate; d)-e) mixed condensate. 
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FIG. 3: E mass with current I. m^ = 0.10 GeV, t = 1 and A = 0.44 GeV. 
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FIG. 4: Leading terms of the R.H.S of (HH) with current I. m^ = 0.10 GeV, t = 1 and A = 0.44 
GeV. Solid line: perturbative term; dotted line: operators of dimension 4; dashed line:operators of 
dimension 6. 
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FIG. 5: Relative strength of the pole (solid line) and the continuum (dotted line) as a function of 
the Borel mass squared with current I. rag = 0.10 GeV, t = 1 and A = 0.44 GeV. 
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FIG. 6: E mass with current II. m, = 0.10 GeV and A = 0.24 GeV. 
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FIG. 7: Leading terms of the R.H.S of (|14j) with current II. 



m.o 



0.10 GeV and A = 0.24 GeV. Solid 



line: perturbative term; dotted line: operators of dimension 4; dashed Iine:operators of dimension 
6. 
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FIG. 8: Relative strength of the pole (solid line) and the continuum (dotted line) as a function of 
the Borel mass squared with current II. nis = 0.10 GeV, t = 1 and A = 0.24 GeV. 
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